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A search for an exotic natural radioactivity of lead nuclei, using nuclear emul- 
sion sheets as detector, is described. We discuss the experimental set-up of 
a test performed at the Gran Sasso National Laboratory (Italy), the event 
simulation, data analysis and preliminary results. 
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1. Introduction 

In the late 80's, some theoretical work was dedicated to investigate possible 
exotic types of nuclear radioactivity, consisting in the emission of light 
particles such as pions or muons from heavy nuclei. 1 ' 2 

Muons or pions could be emitted by nuclei through the decays: 2 

(A,Z) ^ n ± + v li {V ll ) + (Ai,Zi) + ... + (A n ,Z n ) (1) 

(A,Z) -^ir ± (ir°) + (A 1 ,Z 1 ) + ... + (A n ,Z n ) (2) 

where, for reasons of energy and momentum conservation, the number of 
fragments n is > 2. The nuclei (Ai,Zi), (A 2 , Z 2 ), (A n , Z n ) would yield 
a sequence of j3~ decays leading finally to stable nuclei with a balanced 
neutron-proton ratio. 

In ref. 2 some nuclear charge thresholds for different possible sponta- 
neous particle emission were listed: 

(i) [i (prompt muon) for Z > 72 

(ii) 7T^~ (prompt pion — > delayed muon) for Z > 76 
(m) 2/x (prompt muon pairs) for Z > 91 
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(iv) 27r ± (prompt pion pairs — > delayed muon pairs) for Z > 100 

Natural lead is mainly composed by three nuclides: 206 Pb (24.1%), 207 Pb 
(22.1%) and 208 Pb (52.4%). They are stable nuclides, but the channels (i) 
and (ii) are energetically allowed. 

In the hypothesis of the decays (1) and (2), the fission fragments would 
remain nearly at rest; most of the available energy would be used to pro- 
duce the fi (or 7r) and the kinetic energies of /x and f M (or it). The total 
kinetic energy in a decay in 2 fragments with close values of A\ and 
A 2 (symmetric fission) is ~30 McV for negative muons and ^20 MeV for 
positive muons. Considering that the associated muon neutrino takes away 
a sizable fraction of this energy, the spectrum of emitted muons could be 
like in a (3 decay, with an average energy of 10 15 McV. 

The spontaneous or neutron induced fission of Pb has never been ob- 
served. A search for a Pb muonic decay can be made as a byproduct of 
the OPERA experiment, 3 aimed to confirm neutrino oscillations in the pa- 
rameter region indicated by some atmospheric neutrino experiments 4 -. 6 In 
the following, we discuss the experimental set-up, a Monte Carlo (MC) 
simulation and the reachable limits. 

2. Experimental set-up 

We propose to perform an experimental search for muonic radioactivity 
from lead nuclei in the low background conditions offered by the Gran 
Sasso underground Laboratory (LNGS). The low cosmic muon flux and 
the low natural radioactivity of the rock in the experimental halls of the 
LNGS provide unique conditions, allowing a potential discovery, or, at least, 
to establish a good upper limit for this exotic decay process. A detailed 
description of the different background sources is given in ref. 7. 

The OPERA base element ("brick") is composed of 56 lead sheets (1 
mm thick) interleaved with 2 nuclear emulsion sheets (43 /im thick) on both 
sides of a 200 /zm thick plastic base. The area of each sheet is 10.3 x 12.8 
cm 2 . The OPERA bricks (each containing a mass of 8.23 kg of Pb) could 
allow an experimental search for muon emission from lead, with exposures 
of several months. Their analyses with the fast automated optical micro- 
scopes 8 would establish the local background contributions and validate the 
analysis procedures. 

As the background rejection/reduction (see ref. 7) is a crucial point for 
this search, the detectors are surrounded on all sides by a shield, making a 
closed box structure. The shield is composed of an inner layer 5 cm thick of 
very pure copper followed by 15 cm of very low activity lead. The third layer 
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of the shield is a 3 cm thick polyethylene, in order to absorb neutrons a . The 
set-up is located in the emulsion storage room, in hall B of the Gran Sasso 
Laboratory. The radon reduction is obtained with a ventilation system with 
fresh air forced circulation. 

3. Monte Carlo simulation 

A MC simulation program was implemented to estimate the occurrence of 
different event topologies. The simulation is based on the GEANT 9 package 
applied to the OPERA lead/emulsion set-up. The simulation reproduces 
one complete OPERA brick, where muons of different energies (see the first 
column of Table 1) originate in random positions in the lead sheets. The 
initial muon directions are isotropically generated. We assumed different 
definitions for a candidate event, requiring that the muon crosses at least: 
(i) 2 emulsion films, (it) 4 emulsion film (and thus also a lead plate) and 
(iii) 6 emulsion films (and thus also 2 lead plates). We also requested the 
detection of the decay positron or electron in at least (iv) 2, (v) 4 or (vi) 6 
emulsion layers, together with the muon detection. 




Fig. 1. Simulated events of spontaneous fi + (left) and fi (right) emission from lead 
inside an OPERA brick, assuming an initial muon kinetic energy of 15 MeV. 

For each topology, the values quoted in Table 1 were obtained as av- 
erages over all possible muon emission directions. These estimates may be 
considered as geometrical efficiencies e g . The percentages of events listed in 
Table 1 were computed for samples of 1500 MC events with fixed energies. 

In Fig. 1 are shown two simulated muon emission events from the lead 
inside an OPERA brick, assuming an initial kinetic energy of 15 MeV, and 
the production of a fi + (left) and a fj,~ (right). 



a Thcse thicknesses have been found to be adequate with a Monte Carlo simulation 
considering the effects of 0.5-^2.6 MeV photons, which could produce electrons mimicking 
the searched events. 



January 21, 2010 16:30 



WSPC - Proceedings Trim Size: 9in x 6in Giorgini'Proc'ICATPP09 



4 



Table 1. Percentage of events as a function of the minimum number of emulsion 
films crossed by both fi (iV M ) and e (N e ), for different fi initial energies (col. 1). 
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4. Estimates of global detection efficiencies 

The geometrical e g have to be multiplied by the and e ± reconstruction 
efficiencies to obtain the total efficiency. With the OPERA tracking proce- 
dure, 10 the mean detection efficiency for each "microtrack" in one emulsion 
film is ~ 95%. The instrumental limit of 0.8 rad on the incident direction 
introduces an event selection of ~30%. The "base track" (obtained from 
2 microtracks separated by the plastic base) reconstruction efficiency is 
~ 90% (95% x 95%). The top-bottom linking efficiency, mainly due to the 
multiple scattering in the plastic base, ranges from the ~50% for 5 MeV 
particles to ^99% for particles with energies > 20 MeV. The base track 
linking efficiency is ^6% for the whole range of initial muon energies. 

The detection efficiency e to t for a muon crossing at least n emulsion 
films is the product of the percentage of events with these topological re- 
quirements e g (col. 6 of Table 1), the 30% given by the angular limit, n 
times the 95% efficiency for the microtracks, n/2 the top-bottom linking 
efficiency and the linking efficiency between the base tracks. Table 2 shows 
the % efficiencies e g and e to t computed requiring at least 2 emulsion films 
crossed by the muon (col. 2-3) and 4 emulsion films crossed by the muon 
(col. 4-5). The values refer to the present OPERA tracking. 10 

As the expected halflifes ii/ 2 are much larger than any reasonable ex- 
posure time T, the expected sensitivities are estimated from 

W = 7 — Tetot (3) 

where SN — 2.3 is the number of events corresponding to a 90% C.L. limit 
assuming no candidates, Nq is the initial number of nuclei, and etot is the 
experimental efficiency. Using one OPERA brick for one year exposure and 
a global detection efficiency of ~10%, a sensitivity of ^7 • 10 23 yr (90% 
C.L.) could be reached. 
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Table 2. Geometric efficiencies (e 9 ) and total muon detection efficiencies 
(ttot) as a function of the initial muon energy computed on the basis of the 
OPERA tracking procedure for = 2 (col. 2-3) and for = 4 (col. 4-5). 
All values are in %. 
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5. Conclusions and perspectives 

A test search for spontaneous emission of muons from Pb nuclei, using 
some OPERA lead/emulsion bricks, was described. We are in the process 
of making a complete simulation of the detector including its response and 
the track reconstruction efficiencies. We have shown that stringent limits 
for spontaneous muon radioactivity may be reached: t\ji > 7 • 10 23 years. 
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